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SUMMARY 


The steady-state flow of an infinitely conducting, rarefied, neutral 
plasma past a line current is considered. According to the Chapman -Ferraro 
theory, the plasma Is confined to the exterior and the magnetic field to 
the Interior of a cayity, the boundary of which is a thin current sheath. 
Within these restrictions the problem is solved exactly by a method similar 
to the hodograph method of Incompressible fluid dynamics, and the magnetic 
field inside the cavity is presented. The boundary so obtained agrees with 
that previously given by Hurley who used a different method. Also, the 
magnetic field inside the cavity is determined by means of an approximate 
solution for the cavity shape previously obtained by Ferraro. Results of 
the approximate and exact solutions are compared* 


INTRODUCTION 


Chapman and Ferraro, in studies of the connection between solar 
activity and geomagnetic storms, have considered the idealized problem of 
a rarefied neutral plasma stream interacting with a permanent magnetic 
field. Under these conditions the permanent field does not penetrate the 
plasma but has the effect of producing a hollow or cavity In the plasma. 

The steady-state formulation of this problem is summarized by Dungey 
(ref. l) . Several Investigators have applied this formulation to the 
interaction of a plasma with a two-dimensional dipole (refs. 2 , 3 , k ) . 

The recent evidence that a ring current flows in the magnetosphere 
has given impetus to the study of the interaction of a plasma with a two- 
dimensional line current (see ref. 5 ). Ferraro (ref. 6 ) and Hurley (ref. 7 ) 
have initiated such studies recently. The former obtained an approximate 
expression for the boundary between the stream and the magnetic field, and 
the latter an exact one. The present paper applies the formulation of 
Hurley and Ferraro to the exact and approximate magnetic fields inside the 
cavities . Inasmuch as the approximate formulation of Ferraro has been 
applied to three-dimensional problems by Beard (ref. 8 ), and Spreiter and 
Briggs (ref. 9 ), extensive comparisons of the exact and approximate 
solutions are made. 
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FORMULATION OF THE PROBLEM 


It Is considered (ref. l) that the particles in the stream travel up 
to the boundary of the hollow in straight lines and are in effect reflected 
specularly. In so doing, each particle imparts to the surface twice its 
normal component of momentum. This change in momentum is balanced by the 
magnetic pressure at the boundary surface. One thus obtains 

= p D cos 2 a (la) 

8rt 


p D = 2mnv2 


(lb) 


where H s is the total magnetic field at the boundary surface; m, n, and 
v are the mass, number density, and velocity of the positive ions in the 
undisturbed stream; and a Is the angle between the stream and the normal 
to the surface as shown in figure 1. Because no part of the surface may 
be shielded from the stream, cos a must lie between zero and unity. 

The condition that the magnetic field does not penetrate the plasma 
means that at the boundary surface the magnetic field must be tangent to 
the surface, that is, 


& = 5E (2) 

dx H x 


where x and y are two-dimensional rectangular Cartesian coordinates 
(see fig. l) and H x and Hy are the x and y coordinates of the magnetic 
field. 

The equations to be satisfied inside the cavity are: 

8H X (jHy , . 

r-* + r-E = 0 (3) 

ox oy 

and. 

dHy BE X 

dx dy 


w 
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The field components must approach the appropriate values at the line 
current, that is. 


H x - 

.21 Z 

(5a) 

r 2 

By - 

2lx 

r 2 

(5b) 


where r 2 = x 2 + y 2 and I is the current carried by the wire. In the 
above the current is considered to flow out of the paper, and is measured 
in abamps, and the field is measured in gauss. 

THE EXACT SOLUTION 


Predictions as to the nature of the physical problem can be made once 
equations (3) and ( 4) have been solved subject to the boundary conditions 
of equations (l), (2 ), and (5) (with the restriction that cos a is 
greater than zero) . The main difficulty in obtaining a solution to the 
equations in their present form is that conditions (l) and (2) must be 
satisfied on an unknown boundary. To state this another way, the boundary 
must be determined so that (l) and (2) are both satisfied on it. It will 
now be shown that, when independent and dependent variables are inter- 
changed, the boundary becomes known. Before this is done, however, it is 
convenient to introduce a potential function. 

If we Introduce the complex function W, that Is, 

W = « + it 

where 0 is a potential function and t a stream function, and define 

H x = — = (6a) 

dx Sy 

H y . 3* = . H (6b) 

6y 6x 

equations (3) and (4) become 

d 2 3> + S 2 $ 

6x 2 3y 2 


(7a) 



4 


(7b) 


dfi + sfi = o 

8x 2 dy 2 

Equations (6) are recognized to be the Cauchy -Riemann equations; thus W 
is an analytic function of the complex variable z = x + iy (see ref. 10 ). 
Similarly, equations (3) and (4) are the Cauchy -Riemann equations for the 
complex function H* = H x - iHy, and H* is also an analytic function 
of the complex variable z. Equations (7) can be combined into 


d 2 W + dfw _ 
8x 2 c3y 2 


(8) 


A result of complex function theory which can be shown by direct substi- 
tution (making use of eqs . ( 3 )? (4) and (6)) is that if x and y are 
replaced by H x and Hy, equation (8) becomes 


8 2 W + 8 2 W = 0 

8H X 2 8H y 2 


(9) 


The boundary conditions, equations (l) and (2), will now be combined into 
a single relation involving only the component of the magnetic field. 

From figure 1 


-cot a = ^ ^ 

dx H x 


Alternatively, 


cos a = 


3l 


JHx 2 + Hy 2 

Substituting into equations (l) and simplifying, we have 


( 10 ) 


( 11 ) 


E x 2 + Hy 2 = ±(8xp o ) l/2 Hy 


(12) 


where the positive sign is chosen because Hy is always positive on the 
boundary. Equation (12) then becomes 


H x 2 + [Hy - (2*P 0 ) 1/2 ] 2 = 2rtp 0 


(13) 
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This is the equation of a circle as shown in figure 2 . The condition 
stated in equation (5) can he expressed in terms of W and the magnetic 
field as follows. Near the wire 


W -*> -2il log z 


(i*0 


where I is the current carried by the wire. Equation ( 5 ) i n complex 
notation is 


( 15 ) 


Substituting into (l 4 ) we have, for large H* 


W - 2il log =S_ (l6) 


The boundary value problem has now been transferred into a fixed boundary 
value problem. It is now necessary to seek a solution to equation ( 9 ) 
which will satisfy the condition that the field is tangent to the curve 
given by equation ( 13 ) an( l 'will satisfy (16) for large H*. It will be 
recognized that this is analogous to a circular flow which at infinity 
approaches a vortex flow. 

It seems desirable , before obtaining this solution, to point out 
aspects of the method which apply to other two-dimensional problems . For 
any enclosed current system, such as a dipole, a relation, such as equa- 
tion (16), can be obtained relating the magnetic field to the complex 
potential W. Equation (13) applies to any two-dimensional problem to 
which (l) and (2) apply. In more complicated problems, however, it Is 
not possible to choose a single sign in this equation as was done above . 
The boundary therefore becomes a circle plus arcs of another circle. Thus 
any two-dimensional problem can be reduced to a fixed boundary value prob 
lem In the "magnetic" plane. It may happen, however, that the mapping 
from the physical plane to the magnetic plane Is not one to one; that is, 
two or more points in the physical plane may map into the same point in 
the magnetic plane. 

The solution to the present problem can now be written by analogy 
with vortex flow, that is, 


1H* - (2itp n ) 1/2 
21 


W = 21 i log 


( 17 ) 
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Since 


— = H* 
dz 


(18) 


an expression Involving only W and z can be obtained by solving equa- 
tion (17) lor H*, equating this to — , and integrating. The result is 

dz 


z = 


hi 


(8»pJ 1/2 


log 


1 + 


(8»pJ 1/2 i £r 


4l 


(19a) 


or, separated into real and imaginary parts, 


x + iy = 


41 


(8itp 0 ) 


1/2 


log 


1 + 


(8flp n ) l/2 - 


COS 


( 8jtPn) 1 

41 


/ 2 _ 4 


4l 

2 1 


21 


21 


sin 


21 


+ i 


41 


( 8?tp 0 ) 1/2 ^ 21 0 


(8*p n ) 


1/2 


tan -1 


41 


sin -r- 

21 


(ea Pg) 1/a e - 21 cos A 


1 + 


41 


21 


(19b) 


Since lines of constant ^ are field lines, they can be determined if 
particular values of ^ are chosen and z Is calculated for various 
values of <P. In particular, the open field line, vhich is the boundary, 
can be determined from 


(8jtPo) 1/2 

4l 



(20a) 


since this produces a curve that approaches infinity in the negative x 
direction. Substitution of (20a) into (19) leads to 
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- j. (x + iy) = 8 log ^cos + i jjj ( 20 b) 

This is identical to the boundary previously determined by Hurley (ref. 7). 
The value of the field function at the nose of the boundary (and indeed on 
the complete boundary) is from equation (20a) 


tyjq- — -21 log 


4l 

(8*p 0 ) 1/2 


(21a) 


The difference in the values of i | r associated with two field lines 
represents the magnetic flux passing between them. Therefore the flux 
passing between any field line and the boundary is 


t =M |r - % 


(21b) 


An expression for the components of the magnetic field can be obtained by 
inverting equation (17) • The result is 


H 


l x 


+ iH v = 2(8:rtp 0 ) 1 / 2 


-sin £ 


(8jtp 0 ) 


41 


x/a _ 
e 


* 

21 


+ i 


cos 


21 


(8mp 0 ) 


41 


x/a _ 
e 


* 

21 


+ 1 


1 ( 22 ) 


The field lines are plotted in figure 3. The variation of the absolute 
value of the magnetic field is indicated by means of ticks. The direction 
of the field is, of course, tangent to the field line. Coordinates of the 
field lines and components of the magnetic field are given in table I. 

It was convenient to use stretched polar coordinates, that is, 

r = (§f§J /2 (* 2 + y 2 ) 1/2 (23a) 


6 = tan" 1 ^ (23b) 

X 


The components of the magnetic field are expressed as: 


Hp = H x cos G + Hy sin 0 
H q = -H x sin 6 + Hy cos 6 


(24a) 

(24b) 



8 


THE APPROXIMATE SOLUTION 


If the equation of the boundary field line is known, the problem of 
determining the magnetic field Inside the hollow becomes a fixed boundary- 
value problem; that is, the boundary conditions can be applied on a curve 
known beforehand. If an approximate boundary is known, it can be used to 
obtain an approximate solution to the problem. It Is reasonable to suppose 
that if the approximate boundary is close to the actual boundary, the 
approximate field' lines will be reasonably close to the actual field lines. 

Ferraro (ref. 7) obtained an approximate boundary by assuming that 
the total field at the boundary is equal to 2f times the component of 
the magnetic field of the wire parallel to the boundary. He then used 
equation (l) to relate the known magnetic field to the slope of the 
boundary. This led to a differential equation for the boundary shape. 
Ferraro's solution Is 


where 


T) « 


cot q 


(25a) 

X 

r o 


( 25 b) 

X 

r o 


(25c) 

2l 2 f 2 


(25<1) 

Ferraro, In his 
given in (lb) . 
then determined 

original development, 
For a discussion of 
f by requiring that 


r o 2 


and p Q is given by equation (lb), 
set p Q equal to one half the value 
this point see reference 9* Ferrarc 
the magnetic field just outside the nose of the hollow be zero. (Since 
the boundary is approximate, the field everywhere outside the hollow cannot 
be set to zero.) He found that f = 0.68 approximately. It Is convenient 
in the present case to determine f so that the approximate boundary 
matches the exact one at the nose. This gives 


f = 0.6931 


(26) 


We now regard the boundary as known and proceed to calculate the 
field inside it. The differential equation for the field lines is 


d Z = 
dx H x 


(27) 


If the components of the magnetic field are knovn, this equation can be 
integrated to determine the coordinates of the field lines. The components 
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can be deter min ed In tvo ways since two boundary conditions exist for the 
exact problem. The first is equation (l) which specifies the value of the 
field on the boundary. The other is equation (2) which specifies that the 
normal component of the field vanishes on the boundary. Associated with 
the given approximate boundary and each of these conditions is a unique 
set of field lines. Therefore two different sets of field lines can be 
computed inside the approximate hollow. As mentioned earlier, Ferraro 
used equation (l) to determine his approximate boundary. It is therefore 
more consistent to use equation (l) in conjunction with Ferraro's boundary 
to determine the field lines. The magnetic field Is made up of the field 
caused by the current flowing in the wire and that caused by the currents 
flawing In the bounding surface. The currents flowing in the sheet can 
be deter min ed by applying equation (l) . Since the magnetic field outside 
the sheet Is zero (nearly zero In the approximate case) and Is equal to 

(8jp 0 ) l/2 cos a just inside the sheet, the currents flowing must be such 
as to produce a jump of this magnitude In the magnetic field, that is, 

. H s /Po If / O o\ 

j = — = - / — cos a = - cos a (28) 

4jc v 2rt *r 0 

where j Is the current per unit length of the sheet. The Biot-Savart 
law can now be used to obtain a relation for the magnetic field induced 
by these currents. The result is 


H 


Xi 


r * r o , 

f dy’ 

Jtt'o J-xr 0 (x - x’) 2 + (y - y') 2 


(29a) 



2lf 


j.„ 0 (x - X ’) 2 + (y - y ') 2 7 


(29b) 


since ds cos a = dy, where ds is an element of arc length along the 
sheet. The total field can be obtained by adding the field of the wire 
to equations (28). If equations (25) are used the field components are: 


Hx = 1 
(Sjpo) 1 / 2 4 f 


ign 

I 2 + q 2 


A* 

2f / 

J-rt (l sin q 


(■q - q ')sin 2 T) ' 


q'cos T)’) 2 + (q - q') 2 sin 2 q' 


dq ' 


(30a) 
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H- 


Z = _L 

(8^p 0 )Va kt 


g~n 
I 2 + fi 2 


2f 


r (| sin tj 1 - tj 'cos tj * ) sin tj 1 


rt 'J-n (l sin tj' - tj 1 cos tj ' ) 2 + (tj — ii , ) 2 sin 2 r| ’ 


dTj ' 


(30b) 


The field lines can be computed by integrating the equation 


dll _ dy _ % 
d£ dx H x 


( 31 ) 


Equation (31) can now be solved by computing H x and Hy from equations (30) 
and integrating numerically. The integrals in equations (30) are free from 
singularities except when the integral curve lies on the bounding sheath. 

In that case the denominators of the integrands vanish. Because of this 
difficulty, the Integrals were split into three parts: 



(tj - Tj * )sln2 , n 1 

(£ sin tj 1 - T) 'cos tj') 2 + (tj 


dt ! 1 

rj ' ) sin 2 Tj ' 



(tj - tj ' )sin 2 Tj ' 

(£ sin tj ' - T| * cos rj') 2 + (tj 


- tj ' ) 2 sin 2 Tj 1 


dr) ' 


r* 

(l sin tj ' 


(tj - tj ' )sin 2 Tj 1 
T) 'cos TJ 1 ) 2 + (tj 


tj ' ) 2 sin 2 rj ' 


dTj ' 


+ Ki 
(32a) 
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(l sin T| 1 - T) 1 cos tj '~)sin tj 1 


(| sin T)' - tj 'cos tj ') 2 + (tj - tj ' ) 2 sin 2 Tj ' 


P (l sin tj 1 - tj ' cos Tj')sin tj' d r 

’-it sin T| 1 - tj'cos tj ’) 2 + (TJ - tj ' ) 2 sin 2 Tj ' 


r 1 (| sin T] ' - tj 'cos Tj')sin tj ' 

T|+€ (i sin tj' - T) ’cos tj') 2 + (n - t) ') 2 sin 2 T) 


- dr] ’ + Ks 
(32b) 


K f (V - abaai an' (32c) 

J-n-e (I sin tj' - tj'cos tj ') 2 + (tj - Tj'^sii^Tj' 


K 2 = / (6 sin n- - t|-coh i|-;bxu h dTJ r 

^rj-e (I sin rj ' - tj 'cos tj ') 2 + (tI - Tj , ) 2 sin 2 Tj' 


‘ T - "'cos n')sin n* 


( 32 d) 


The singularities are now contained in K]_ and K 2 . They were then 
evaluated approximately as follows: The substitution 

Tj' = T| + 5 

was made in equations (32c and (32d), and linear terms in 8 were retained. 
This approximation amounts to replacing the actual boundary in the neigh- 
borhood of tj 1 = tj with a straight segment whose slope is equal to the 
boundary slope at tj. The resulting integrals can then be evaluated 
analytically. Thus Ki and K 2 were found to be approximately 

Kl = 1 log < A - s + Be ) a ± - B tan -i B(A - il 1 e(B2 + 1) 

2(B 2 + 1) (A - | - Be) 2 + e2B2 + l A - | 


JL_ tan -i B(A - j) - e(B ± 1), 


B 2 + 1 


A - i 



12 


Ife * s log * Sill +. A? + 1 tan-l gl A — J) * 12 

2(B 2 + l) (A - i - Be) 2 + £2 B 2 + 1 A - | 


- 1 tan- 1 B(A - I) - e(B° t l) 

B 2 + 1 A - | 


(33b) 


where 


A = 7) cot T| (33c) 

B = cot T] - - \ (33d) 

sin^T] 

Equation ( 31 ) was then integrated numerically. An approximate field line 
was computed corresponding to each exact field line obtained previously. 

In every case except for the field line starting at the nose of the 
boundary, the starting point for the integration was chosen so that the 
exact and approximate curves coincided on the positive x axis. When 
integration was started precisely at the nose of the boundary, numerical 
difficulties arose. It was found, however, that Integration could be 
started just Inside the boundary. The starting value chosen was 0. 99999 
times the x coordinate of the nose . Small variations In the starting 
point did not affect the resulting field line appreciably. Results of 
the calculation are compared with the exact solution in figure k and are 
given in more detail In figure 5* In figure 5 the magnitude of the mag- 
netic field Is indicated by means of ticks. The magnetic field strengths 
Indicated on the bounding field line were calculated just inside the 
boundary. Numerical results for coordinates and magnetic field components 
are given In table II. As In the exact case it was convenient to use 
polar coordinates. They are related to £ and tj by 

r = fU 2 + n 2 ) l/2 (34a) 

0 = tan -1 (34b) 


The components of the magnetic field, Hjr and H q, are related to H x and 
Hy by equations (24) . 

The magnetic flux passing between the nose of the boundary and any 
field line can be computed as follows. From equation (6b) 


x N 


Hy dx * 


f = ^ - % = 


(35) 
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where xjj is the x coordinate of the boundary nose and x is any 
point in the field. Substituting from equation (25b) gives 


t = 


-4lf 


r_st 

(8jtp 0 ) 1/2 


d^» 


(36) 


The integration was carried out numerically for each field line computed 
and the results are given in table II. 


RESULTS AND DISCUSSION 


An exact solution for the interaction of a neutral rarefied plasma 
stream with the magnetic field of a line current has been presented. 
Although the shape of the boundary has been obtained previously, the 
magnetic field is presented here for the first time. The method employed 
is applicable in part to many such two-dimensional problems. The magnetic 
field inside Ferraro T s approximate boundary has also been determined. It 
follows from a comparison of the boundaries that the numerical value for 
f obtained by Ferraro is extremely close to the value one obtains by 
matching the approximate with the exact boundary at the nose. From 
figure 4 it can be seen that exact and approximate field lines which begin 
at the same point remain extremely close together in the upstream portion 
of the figure . In the downstream portion the most serious divergence 
between pairs of field lines is that between the exact bounding field line 
and the approximate field line which almost coincides with it at the nose. 
In this case the approximate field line forms a closed curve while the 
exact one goes to infinity in the negative x direction. Since the 
approximate field line is not altered significantly by slight variations 
in the starting point of -the integration, it is concluded that this diver- 
gence is not due to a numerical difficulty but is Inherent in the approxi- 
mation. It can be seen from figures 3 and 5 (or tables I and II ) that in 
the upstream part of the figure the magnetic field strengths on corre- 
sponding field lines are very nearly the same. The disagreement in the 
downstream part of the figure is due mainly to the divergence of the field 
lines themselves . 

It has been shown above that the approximate field lines, with the 
exception of the one starting at the nose, agree with the exact ones to 
the same order (but opposite sense) that Ferraro f s approximate boundary 
agrees with the exact one. The disagreement of the approximate field 
lines with the exact ones Is a result of the disagreement of the approx- 
imate boundary with the exact one and not a result of the numerical 
confutation scheme . 


Ames Research Center 

National Aeronautics and Space Administration 
Moffett Field, Calif., Aug. 10, 1962 
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TABLE I . - EXACT SOLUTION FOR MAGNETIC FIELD LINES AND FIELD COMPONENTS INSIDE SHEATH 
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TABLE II . - APPROXIMATE SOLUTION FOR MAGNETIC FIELD LINES 
AND FIELD COMPONENTS INSIDE SHEATH 
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Figure 2 .- The sheath in the magnetic plane. 


















